arXivil507.00199v2 [cond-mat.mes-hall] 28 Oct 2015 


Nanoscale phase-engineering of thermal transport with a Josephson heat modulator 


Antonio Fornieri, 1, * Christophe Blanc, 1 Riccardo Bosisio, 2,1 Sophie D’Ambrosio, 1 and Francesco Giazotto 1, ^ 

l NEST, Istituto Nanoscienze-CNR and Scuola Normale Superiore, Piazza S. Silvestro 12, 1-56127 Pisa, Italy 
2 SPIN-CNR, Via Dodecaneso 33, 1-16146 Genova, Italy 
(Dated: October 29, 2015) 


Macroscopic quantum phase coherence has one of its 
pivotal expressions in the Josephson effect [1], which man¬ 
ifests itself both in charge [2] and energy transport [3- 
5], The ability to master the amount of heat transferred 
through two tunnel-coupled superconductors by tuning 
their phase difference is the core of coherent caloritron- 
ics [4-6], and is expected to be a key tool in a number of 
nanoscience fields, including solid state cooling [7], ther¬ 
mal isolation [8, 9], radiation detection [7], quantum in¬ 
formation [10, 11] and thermal logic [12]. Here we show 
the realization of the first balanced Josephson heat mod¬ 
ulator [13] designed to offer full control at the nanoscale 
over the phase-coherent component of thermal currents. 
Our device provides magnetic-flux-dependent tempera¬ 
ture modulations up to 40 mK in amplitude with a maxi¬ 
mum of the flux-to-temperature transfer coefficient reach¬ 
ing 200 mK per flux quantum at a bath temperature of 
25 mK. Foremost, it demonstrates the exact correspon¬ 
dence in the phase-engineering of charge and heat cur¬ 
rents, breaking ground for advanced caloritronic nanode¬ 
vices such as thermal splitters [14], heat pumps [15] and 
time-dependent electronic engines [16-19]. 

When two superconductors (Si and S 2 ) are coupled by 
means of a thin insulating layer (I), they form a Josephson 
junction (JJ) through which a dissipationless charge current 
7j can flow due to Cooper pair tunneling. 7j obeys the well- 
known expression [1] 


7j=7osin<p, (1) 

where 7o is the critical current of the JJ and (p the phase differ¬ 
ence between Si and S 2 . Although the Cooper pair conden¬ 
sate carries no entropy under static conditions and the super¬ 
current 7j cannot contribute to the heat current [3, 20], it was 
shown [3-5] that the Josephson effect has a thermal counter¬ 
part. Indeed, if we impose a temperature bias to the JJ by 
setting the temperature in Si to be T\ > T2 (T2 being the tem¬ 
perature in S 2 ), a finite stationary thermal current 7 t ot will flow 
through the junction [3, 20, 21], 

Jtot(Tl, T2) = Jqp(T\ , 7^) — 7int(Tl 5 12 ) COS<p. (2) 

In Eq. (2), the term 7 qp accounts for the energy transferred 
by quasiparticles whereas J mi represents the phase-dependent 
component of the heat current which originates from energy¬ 
carrying tunneling processes involving transfer of Cooper 
pairs [3, 20] (see Methods for further details). 


* antonio.fornieri@sns.it 
^ francesco.giazotto@sns.it 



FIG. 1. Quantum modulator structure, a. Pseudo-color scan¬ 
ning electron micrograph of the double-loop heat interferometer. 
Source and drain N electrodes, depicted in green, are composed of 
Al 0 . 98 Mn 0.02 and are connected to the A1 island Si (represented in or¬ 
ange) and to a set of A1 wires (dark cyan) serving either as heaters or 
thermometers. Si is tunnel-coupled to S 2 (Al, dark cyan) by means 
of three parallel JJs forming the double-loop SQUID and to a super¬ 
conducting probe S 3 (Al, dark cyan). All the junctions in the struc¬ 
ture are implemented through A10 x tunnel barriers. The areas of the 
loops are of about 15 pm 2 and 30 pm 2 , b. Schematic picture of the 
device. The setup for electrical and thermal measurements are rep¬ 
resented in black and blue, respectively, while the red lines indicate 
the magnetic fluxes and <I >2 piercing the loops. 


Here we experimentally demonstrate that a double-loop 
direct current superconducting quantum interference device 
(DC SQUID) [22] with three parallel JJs [13] represents a fun¬ 
damental building block to achieve full phase-engineering of 
electronic heat currents at the nanoscale. Our device is ro¬ 
bust against structure asymmetries, providing an almost com¬ 
plete annihilation (up to 99%) of J- mi , whereas a conventional 
single-loop DC SQUID would suffer from a reduced contrast 
in the modulation. Yet, it allows to manipulate at will the heat 
current-phase relation, enabling the generation of exotic ther¬ 
mal interference patterns with large modulation amplitude, 
and enhanced sensitivity to magnetic flux variations [13]. 

The double-loop Josephson heat modulator has been fabri¬ 
cated by electron-beam lithography, three-angle shadow-mask 
evaporation of metals and in situ oxidation (see Methods). It 
basically consists of five different parts, as shown in Fig. la. 
Aluminum (with a critical temperature T c ~ 1.4 K) was used 
to implement all the superconducting elements of the struc¬ 
ture, while Alo. 98 Mno .02 was employed as the normal metal 
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(N) [5, 8 , 23]. The source and drain N electrodes act as ther¬ 
mal reservoirs and are tunnel-coupled to a superconducting 
island (Si) through junctions with normal-state resistances 
R s ~ 5.6 kO and R& ~ 6.25 k£2, respectively. Si defines 
the upper branch of the heat modulator and is connected to 
the lower branch (S 2 ) by means of three tunnel junctions 
forming the SQUID characterized by an equivalent resistance 
^squid — 700 O. The ratio between the areas of the interfer¬ 
ometer loops is ^ 2. A superconducting probe (S 3 ) is also 
connected to Si via a tunnel-junction with a normal-state re¬ 
sistance R p ~ 600 f l . Finally, the source and drain electrodes 
are tunnel-coupled to a set of external superconducting wires, 
forming NIS junctions of area ^ 150 x 200 nm 2 and normal- 
state resistance of ~ 20 k Q each, which can be used as Joule 
heaters and thermometers [7] . 

To prove the perfect correspondence in the modulations of 
charge and heat currents, we performed two kinds of measure¬ 
ments, whose schematics are displayed in Fig. lb. First, we in¬ 
vestigated the dissipationless electrical transport across the in¬ 
terferometer via the superconducting electrode S 3 with a con¬ 
ventional four-wire technique. The junction between Si and 
S 3 was specifically designed to have a larger critical current 
than the SQUID one, thereby allowing an exact determination 
of the interferometer critical current. Voltage-current charac¬ 
teristics were recorded for different values of the magnetic- 
flux threading the loops. A well-defined Josephson current 
was observed with a maximum critical current /o,max — 495 
nA at a bath temperature 7b at h = 25 mK. The magnetic-flux 
interference pattern of the SQUID critical current 7 q for se¬ 
lected values of 7b at h is shown in Fig. 2a. Its behavior can be 
understood by writing the expression of the total Josephson 
current flowing through the SQUID and imposing flux quan¬ 
tization constraints for both loops [ 2 ] : 

fj ,squid = 7a sin(<pA) +7 b sin(<pe) ±/c sin(<pc), 

2 7T <T>i ^ 

(pA-<pB--^+2n7T, (3) 

27T<P 2 „ 

<Pc = <Pb + —Z- vlnn. 

<p 0 

Here, I\ and (p\ are the critical current and the phase difference 
in the i-th JJ (being i = A,B,C). d>i and <1 > 2 are the magnetic 
fluxes in the loops, while <t>o — 2 x 10“ 15 Wb is the supercon¬ 
ducting flux quantum. 

For given <t>i and <£ 2 , Iq is obtained by maximizing /j ,squid 
with respect to (p leading to the following expression: 

k = h \Jl + r\ ± r\ ± 2r x a ± 2r 2 /3 ± 2nr 2 y 
= 7b 

where a = cos(27rd>i/<I>o), /3 = cos(27rd> 2 /<t>o), 7 = 
cos[2^(d>i +<J> 2 )/d>o], r\ =/ A // B =R b /R a and r 2 = Iq/Ib = 
Rb/Rc (see Fig. la). 

Figure 2b shows the excellent agreement between the ex¬ 
periment and the model, which allows us to extract pre¬ 
cise values of the interferometer structural parameters. From 
the theoretical fits we obtain d > 2 = (1.93 ±0.02) d>i, r\ — 
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FIG. 2. Electrical response of the interferometer, a. SQUID to¬ 
tal critical current Iq vs magnetic flux <I>i piercing the smallest of 
the interferometer loops for selected values of the bath temperature 
Tbath- b. Zoom of two /oO^i) characteristics measured at 25 mK 
and 1 K. The full circles are experimental values, whereas the solid 
black lines are theoretical results obtained from Eq. (4). c. Normal¬ 
ized critical current /o/7o, max behaviour in a selected range of <I>i for 
Jbath = 25 mK. The maximum observed suppression of Iq is ~ 99 %. 
The inset shows the maximum critical current /o,max = Io(&t — 9) 
(filled symbols) as a function of Tbath together with the theoretical 
Ambegaokar-Baratoff expectation (red dashed line, see text for de¬ 
tails). 


(0.93 ±0.04) and r 2 = (0.70±0.03). These values are con¬ 
sistent with the geometrical parameters expected from scan¬ 
ning electron micrograph of the device (see Fig. 1). The 
temperature-dependence of /o,max is plotted in the inset of 
Fig. 2c along with the Ambegaokar-Baratoff prediction [24] 
for a JJ with a normal-state resistance J^squid and T c ~ 1.4 K. 
The good agreement between theory and experiment confirms 
the ideal behaviour of our double-loop Josephson heat modu¬ 
lator. Moreover, as shown in Fig. 2c, an optimal suppression 
of the critical current of ~ 99 % is obtained for |<E>i | ~ 10 <t>o, 
highlighting the excellent phase control of the Josephson cur¬ 
rent offered by our device. 

On the thermal side, the interferometric properties of our 
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FIG. 3. Thermal response of the Josephson heat modulator at 25 mK. a. Magnetic flux modulation of the drain temperature /drain for 
several values of the source temperature 7h 0t . b. Thermal model outlining the relevant heat exchange mechanisms present in our interferometer. 
Arrows indicate heat current directions for the operating device, i.e., when T^ ot > T\ > T^ m > /b at h (see text). The only existing magnetic- 
flux-dependent thermal current is /mt, squid (labelled in red), c. Experimental T^ram pattern (filled circles) around 4>i = 0 for two selected 
values of /h ot along with the theoretical results from the thermal model (solid black lines), d. Flux-to-temperature transfer function T versus 
<f>l for the same values of T^ot as in panel c. 


heat modulator were investigated by imposing a temperature 
gradient across the device. We emphasize that our experi¬ 
ment is focused on the heat carried by electrons only. We 
assume that lattice phonons present in every part of our struc¬ 
ture are fully thermalized with the substrate phonons resid¬ 
ing at 7b at h, thanks to the vanishingly small Kapitza resistance 
between thin metallic films and the substrate at low tempera¬ 
tures [4, 5, 8 , 25]. When we inject a Joule power (/i n ) through 
a couple of superconducting probes connected to the source 
electrode (see Fig. la), we can raise its electronic temperature 
Thot significantly above Tbath [25]. The quasiparticle tempera¬ 
ture 7j in Si is therefore increased, yielding a thermal gradient 
across the SQUID. This hypothesis is expected to hold be¬ 
cause the lower branch of the interferometer (S 2 ) extends into 
a large-volume lead, providing efficient thermalization of its 
quasiparticles at /b at h. This thermal gradient gives rise to the 
appearance of a finite heat current / t0 t flowing through each JJ 
of the SQUID, leading to a periodic magnetic modulation of 
the drain electronic temperature (/drain)- Both 7h 0 t and /d ra i n 
are monitored by measuring the voltage drop across two pairs 
of current-biased NIS junctions [7], as shown in Fig. lb. 

Figure 3a shows 7d ra i n oscillations as a function of the mag¬ 


netic flux for different values of 7h ot at 25 mK. Even at a 
glance, the thermal magnetic interference pattern perfectly 
mimics the Josephson critical current one. As 7h 0t becomes 
higher, the average value of /d ra i n ((Zd ra in)) raises well above 
Tbath due to the increased heat flow across the structure, while 
the range of /d ra i n spanned by the oscillations decreases from a 
maximum of ~ 40 mK and tends towards saturation for larger 
Thot- We stress that so large temperature modulation ampli¬ 
tudes have never been reported so far, and stem from the de¬ 
sign choice and implementation of the Josephson heat modu¬ 
lator. 

In order to account for our observations, we have formu¬ 
lated a thermal model outlining all the relevant heat exchange 
mechanisms present in the structure. The model is sketched in 
Fig. 3b, where the terms / source and /drain denote the quasipar¬ 
ticle heat currents exchanged by Si with the source and drain 
electrodes, respectively. The upper branch of the SQUID, as 
previously mentioned, can release energy to S 3 and S 2 by 
means of / pro be an d /squid [see Eq. (2)], but only the latter 
term can be phase-controlled thanks to the double-loop geom¬ 
etry. Finally, electrons in the source and drain electrodes can 
exchange energy at power / e - p h with lattice phonons residing 
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FIG. 4. Behaviour of the device at different bath temperatures, a. Magnetic-flux dependent T^-ain modulations for different values of 
Tbath at ^hot — 700 mK. b. Measured T^ram (full circles) vs. compared to the theoretical fit (solid black line) for T^ot ~ 700 mK and 
7b £l th = 200 mK. c Average drain temperature (7^ ra j n ), maximum temperature modulation amplitude A7drain = — 7^™ , and maximum 

flux-to temperature transfer coefficient |T| max vs. T^ot for different values of T^ath- 


at 7b a th (see Methods for details). 

The steady-state electronic temperatures T\ and 7d ra i n can 
be calculated for any given 7h ot , 7b at h and magnetic fluxes 
T>i and <l >2 by solving the following system of energy-balance 
equations: 

7source = J probe + «TSQU ID (^1,^2) +/drain > (5) 

/drain = /s—ph, drain • (6) 

These equations account for the detailed thermal budget in 
Si and drain electrodes by equating the sum of all incoming 
and outgoing heat currents. In the model we neglect photon- 
mediated thermal transport owing to poor impedance match¬ 
ing between source and drain electrodes [26-28], as well as 
electron-phonon coupling in Si due to its reduced volume and 
low experimental 7b at h [29] . 

By imposing conservation of the circulating supercur¬ 
rent [13] together with flux quantization [Eqs. (3)] for both 
loops, we get 

/SQUID = (1 + ri — (7) 

where /® p and are detailed in the Methods. From Eq. (7) 

it is evident that the expected magnetic flux dependence of 
the heat current is exactly the same obtained for the Joseph- 
son critical current. The calculated 7d ra i n was fitted to the 


measured data by using the structure parameters determined 
from the electrical measurements, and by adding a single phe¬ 
nomenological prefactor to reduce the amplitude of (see 
Methods for further details). 

As shown in Fig. 3c, the model provides a remarkable 
agreement with the experimental data, grasping the relevant 
heat transport mechanisms in our system. Furthermore, this 
theoretical picture allows us to estimate the suppression of the 
phase-coherent component of the heat current in the Joseph- 
son thermal modulator, otherwise not attainable from any 
direct measurement. The fit reveals a maximum contrast 
Wnf X — ^irrt n )/^int^ — 99%, exactly as in the case of To. 

Figure 3d displays the flux-to-temperature transfer coeffi¬ 
cient 7 = d 7dr a in/ <9<I>i for different values of 7h 0t . As pre¬ 
viously mentioned, the double-loop geometry allows us to 
achieve non-trivial heat current-phase relations [13], leading 
to sharp features in 7 and to sensitivities which can be as large 
as ^ 200 mKAFo- We note that the latter value is larger by 
more than a factor of three than that obtained in previous ex¬ 
periments [4, 5]. 

The heat modulator dependence on the bath temperature 
is summarized in Fig. 4. Figure 4a displays 7d ra i n (<F) at 
Thot — 700 mK for different values of 7i, a th- Even though the 
distinctive shape remains always well recognizable, the am¬ 
plitude of the oscillations is progressively depressed by the in- 
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creasing electron-phonon coupling, which affects mainly the 
N drain electrode. Nevertheless, our thermal model still pro¬ 
vides good theoretical fits for the experimental data, as shown 
in Fig. 4b. Finally, Fig. 4c assembles the overall behavior of 
the heat modulator. Together with (/drain), also the maximum 
temperature modulation amplitude AT^ain = 7^™* — 7 d ™ n n and 
the maximum value of |T| are displayed for several values of 
7h 0 t and 7b at h- As previously noted, for low values of 7b at h rais¬ 
ing 7h 0 t produces an initial increase of (/drain) joined to a re¬ 
duction of A7d ra i n and |T| max . All these quantities tend towards 
saturation as 7h 0 t exceeds ~ 600 mK due to the impact of the 
electron-phonon coupling. This effect is even more noxious 
as /bath is increased, but not sufficiently strong to completely 
suppress the magnetic oscillations. The latter are indeed de¬ 
tectable up to 7b a th — 400 mK. 

In summary, we have realized a double-loop Joseph- 
son thermal modulator that provides a complete phase¬ 
engineering of the electronic heat current at the nanoscale. 
Despite the unavoidable presence of asymmetries in the JJs 
composing the SQUID, our device offers a suppression of 
the phase-coherent component of the heat current up to 99%, 
leading to a maximum swing in the temperature oscillations 
of 40 mK and to a flux-to-temperature transfer coefficient 
as large as 200 mKM>o at 25 mK. Yet, the interferometer 
demonstrates the perfect correspondence in the manipulation 
of charge and heat transport, paving the avenue to the ideation 
of exotic coherent caloritronic nanodevices where thermal 
currents can be engineered at will. In this perspective, our 
heat modulator could be furnished with two independent su¬ 
perconducting on-chip coils that would provide separate con¬ 
trol knobs for the magnetic fluxes. This would allow to mas¬ 
ter independently the phase-biasing in the loops [13], open¬ 
ing the way to perform closed cycles in the parameter space, 
and leading to the realization of heat pumps [15] and time- 
dependent quantum heat engines [17-19]. This approach can 
be extremely efficient thanks to the accuracy and high fre¬ 
quency (up to several GHz) at which magnetic fluxes can be 
driven in the system [16, 18]. On the other hand, our quantum 
device might represent the building block at the core of ther¬ 
mal splitters [14] able to control the amount of energy trans¬ 
ferred among several terminals residing at different temper¬ 
atures. Finally, further studies are needed to extend present 
results to higher temperatures, exploiting, for instance, high- 
T c superconductors. 


METHODS 

The devices were fabricated with electron-beam lithogra¬ 
phy and three-angle shadow-mask evaporation of metals onto 
an oxidized Si wafer through a bilayer resist mask. The evap¬ 
orations and oxidation were made using an ultra-high vacuum 
electron-beam evaporator, which allowed us to deposit first 
15 nm of Alo. 98 Mno .02 at an angle of 40° to form the source 
and drain N electrodes. Then the sample was exposed to 800 
mTorr of O 2 for 5 minutes to realize the thin insulating layer 
of A10 x forming the tunnel-barriers in all the SIN junctions 
present in the system. Afterwards, 20 nm of A1 were evapo¬ 


rated at 0° to deposit the superconducting island Si, together 
with the probes used as heaters and thermometers. Another 
exposition of 200 mTorr of O 2 for 5 minutes was required 
to realize the insulating layers of the SIS junctions forming 
the double loop SQUID and connecting the probe S 3 to Si. 
Finally the sample was tilted at an angle of 30° and a depo¬ 
sition of 40 nm of A1 was performed to implement the lower 
branch of the SQUID S 2 and the probe S 3 . The source and 
drain electrodes have a volume V source = 1.1 x 10“ 19 m -3 and 
^drain = 7.5 x 10 _2 ° m -3 , respectively, whereas the volume of 
Si is V Sl = 2 x 10 -19 m -3 . 

All the measurements have been performed in a filtered 
dilution refrigerator down to 25 mK. The thermometers 
were current biased through battery-powered floating sources, 
while the heaters were piloted with voltage biasing in the 
range 0-3 mV, corresponding to a maximum power of ^ 120 
pW injected in the source electrode. Thermometer bias cur¬ 
rents were varied from 5 pA to 100 pA in order to maximize 
the sensitivity in different ranges of temperature and to re¬ 
duce the effects of self-heating and self-cooling [7]. Voltage 
and current were measured with standard room-temperature 
preamplifiers. 

In the thermal model, /source(/hot,Zi) = /Nis(Zhot,Zi) and 

/drain (Zl, 7drain) = /nIS (Zl >/drain), where /nIs(/n,Zs) = 

(2 /e 2 R x ) /°°£K(£,7s)[/(£,7n) - /( e ,r s )]d£, being 

Jo 

f(£,T) = [1 + exp(£/kBZ )] -1 the Fermi-Dirac distribu¬ 
tion, N(e,r) = |9t[(e + /T)/ v /(e + /T) 2 -A 2 (r)]| the 
smeared (if T ^ 0) normalized Bardeen-Cooper-Schrieffer 
density of states (BCS DOS) in the superconductor [30], 
A(T) the temperature-dependent energy gap [2], R x the 
tunnel junction normal-state resistance, e the electron 
charge and k b the Boltzmann constant. On the other hand, 
in the expression of /squid (see Eq. 7) two components 
appear: the first accounts for the heat carried by quasipar- 

f*oo 

tides, /qp(7i,7i, ath ) = (: 2/e 2 R x ) / £K 1 (£,7’ 1 )K 2 (£,r bath ) 

J 0 

[/(£, 7j) —/(£,7b a th)]d£, where Ni and N 2 are the normal¬ 
ized BCS DOS of the superconductors forming the Josephson 
junctions. The second component, instead, represents 
the phase-coherent part of the heat current [3, 20], which 

reads .J bdtb ) = ( 2 /e 2 R x ) J eM l (e,T l )M 2 (£,T bat h) 

[f{e,T 1 ) - / (e,?bath)](te. Here, M n (£,r) = 

|3 [—iA n (T)/y /(e + iT n ) 2 — A 2 (T)]| is the Cooper pair BCS 
DOS in the n-th superconductor, with n = 1,2. In the same 
way, / probe ( Zi, /bath ) = /qp ( Zl, /bath ) — /int ( Zi, /bath ) ■ Finally, 
the electron-phonon coupling in the drain electrode generates 
the term /e—ph, drain (/drain? /bath) = ^^drain (/drain — /bath)’ 

where £ = 4.5 x 10 9 WK“ 6 m “ 3 is the material-dependent 
electron-phonon coupling constant [5, 8 , 23]. 

The theoretical curves in Figs. 3 and 4 have been calcu¬ 
lated from Eqs. (5, 6 ), using the the measured values of R s , 
Rd, R p , quid , n, ? 2 , p = d> 2 /d>i and T^5x 10 _ 3 A(0) 
as determined from the electrical characterization of the de¬ 
vices. They were fitted to the experimental data by introduc¬ 
ing a single phenomenological factor A ~ 0.6 reducing the 
coherent component of the heat current in Eq. 7. Its origin 




6 


may be ascribed to fluctuations of the chemical potential in the 
superconducting island Si, which can induce time-dependent 
phase evolutions. Indeed, in the thermal measurement setup 
(see Fig. lb), our device seems more sensitive to such fluc¬ 
tuations because Si is connected to floating leads S 2 and S 3 
through tunnel barriers. By contrast, in the electrical con¬ 
figuration, we set the phase difference of each JJ present in 
the device by current-biasing our interferometer. In the latter 
case, fluctuations turns out to be negligible, as demonstrated 
by the good agreement between the experimental data and the 
Ambegoakar-Baratoff prediction (see the inset of Fig. 2c). 

Furthermore, we estimated numerically the effect of finite 
self and mutual inductances of the loops. These inductances 
renormalize the values of the magnetic fluxes piercing the 
SQUID [2] but their impact on the behaviour of our device 
turns out to be negligible. 


We also checked that the inclusion of a finite Kapitza re¬ 
sistance in the energy-balance equations produces a variation 
of ~ 0.5% in the calculated drain temperature. Similarly, the 
heat exchange due to the quasiparticle-phonon coupling [29] 
in Si was also neglected, since this contribution results to be 
at least ten orders of magnitude smaller than the thermal flux 
flowing from the source to the drain in the explored range of 
7h 0 t and 7b a th- Finally, the power generated by the electron- 
photon coupling between the source and the drain [27] has 
been estimated to be four orders of magnitude smaller than 
^source and /drain> thereby leading to insignificant corrections 
of our results. 
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